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1 Introduction 



It is well known that the Galilei group G(l, 3) and its representations play the same 
role in non-relativistic physics as the Poincare group P(l, 3) and its representations 
do in the relativistic case. In fact the Galilei group and its representations form 
a group-theoretical basis of classical mechanics and electrodynamics. They replace 
the Poincare group and its representations whenever velocities of bodies are much 
smaller then the speed of light in vacuum. On the other hand, structure of subgroups 
of the Galilei group and of its representations are in many respects more complex 
than those of the Poincare group and therefore it is not perhaps so surprising that 
representations of the Poincare group were described in [1] almost 15 years earlier than 
representations of the Galilei group [2] in spite of the fact that the relativity principle 
of classical physics was formulated by Galilei in 1632, i.e., about three centuries prior 
to that of relativistic physics by Einstein. 

An excellent review of representations of the Galilei group was written by Levy- 
Leblond [3]. It appears that the Galilei group, as distinct from the Poincare group, 
has besides ordinary representations also projective ones (see [2] and [4] respectively). 
However its subgroup - the homogeneous Galilei group HG(1, 3) which plays in non- 
relativistic physics the role of the Lorentz group in the relativistic case - has a more 
complex structure so that its finite-dimensional indecomposable representations are 
not classifiable (for details see [5]). And they are the representations which play a 
key role in description of physical systems satisfying the Galilei relativity principle! 

An important class of indecomposable finite-dimensional representations of the 
group HG(1,3) was found and completely classified in [5]. It contains all repre- 
sentations of the homogeneous Galilei group which, when restricted to its rotation 
subgroup, decompose to spin 0, 1/2 and 1 representations. It was explained in [5] 
and [6] how these representations can be obtained from those of the Lorentz group 
by means of the Inonii-Wigner contractions. 

The representations classified in [5] are interesting from both physical and math- 
ematical points of view. Physically, the related spin values exhaust all ones observed 
experimentally for stable particles. Mathematically, the results of paper [5] are unim- 
provable in the sense that the problem of classification of the representations with 
other spin content is unsolvable in general. This gives us undreamt possibilities to 
describe various non-relativistic (quantum-mechanical and field-theoretical) systems 
of interacting particles and fields with spins 0, 1/2 and 1. For instance, it was possible 
to find the most general Pauli interaction of the Galilean spin-1/2 particles with an 
external electromagnetic field, see [5]. 

Starting with indecomposable representations of the group HG(1,3) found in [5] 
an extended class of linear and non-linear equations for Galilean massless fields has 
been found in our previous paper [7]. In particular, Galilean analogues of the Born- 
Infeld and Maxwell-Shern-Simon systems were discussed there. Moreover, all possible 
Galilei-invariant systems of first order equations for massless vector and scalar fields 
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had been classified in [7]. 

In the present paper we continue and complete the research started in [5]- [7]. 
Namely, we study first vector and spinor representations of the homogeneous Galilei 
group in detail and then using them we construct various wave equations for massive 
particles with spin 0, 1/2, 1 and 3/2. In this way a completed list of systems of Galilei 
invariant first order partial differential equations is found which describe particles with 
spin s < 3/2. 

There are well-developed methods for deriving relativistic wave equations invari- 
ant w.r.t. the Poincare group 1 . Consequently several systematic approaches to rela- 
tivistic equations are now available, the most popular of which are those by Bhabha 
and Gelfand and Yaglom [9], those by Bruhat [10] and those based on Garding's 
technigue [11]. These approaches can be used for construction of the Galilei- invariant 
equations as well. 

First, we begin with the Bhabha and Gel'fand- Yaglom [9] approach which is a 
direct extension of the method yielding the Dirac equation. The corresponding rela- 
tivistic wave equations can be written as systems of the first order partial-differential 
linear equations of the form: 

(/3X + £«*(x,0 = 0, (1) 

where p° = ij^, p a = i-£r (a, = 1,2,3), and (3^ (// = 0,1,2,3) and are square 
matrices restricted by the condition of the Poincare invariance. Notice that in the 
relativistic approach matrix (3^ is usually assumed to be proportional to a unit matrix. 

The theory of the Poincare-invariant equations (1) is clearly explained in detail 
for instance in the Gel'fand-Minlos-Shapiro book [12]. 

The other approaches make use of tensor calculus, and the associated equations 
have the form of covariant vectors or tensors (see, e.g., [13]). A popular example is 
the first order Proca equation [14] , 

Vv W^ = re** 4 , ^ ' 

where ^ and is a four-vector and a skew-symmetric spinor respectively which 
transform according to the representation -D(|, \) ® D(l, 0) © D(0, 1) of the Lorentz 
group. The second order Proca [14], Rarita-Schwinger [15] and Singh-Hagen [16] 
equations serve as other examples. Let us mention that all these relativistic equations 
violate causality or predict incorrect values for the gyromagnetic ratio g. The tensor- 
spinorial equations for particles with an arbitrary half-integer spin which are not 
violating causality and admit the right value for g were discussed in detail in [17]. 

1 These methods have had a long a tangled history. It began with Schrodinger's proposal of a 
relativistic equation in 1926, Diracs equation for massive spin 1/2 particle in 1928 and Majorana's 
pioneering study in 1932 and has always been a subject of interest of many scientists (Bhabha, 
Harish-Chandra, Wild, Ficrz, Pauli, Kcmmcr, Duffin, Bargmann, Wigner, Weinberg, Schwingcr, 
and so on, see, e.g., [8] and references cited therein) 
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In the present paper we use both the above mentioned approaches and derive the 
Galilei-invariant equations for particles with spins 0, |, 1 and §. Moreover we present 
a complete list of the Galilei-invariant equations (1) for scalar and vector fields. In 
contains equations obtained via contractions of known relativistic equations as well as 
new Galilei-invariant equations which are received via contractions of new relativistic 
equations. 

A special class of Galilei-invariant equations for particles with arbitrary spins s 
was described in the significant paper [18]. The related equations have the minimal 
number of components and predict the same value g — 1/s for the gyromagnetic ratio 
as it is done by the relativistic wave equations [18]. 

Some particular results associated with the Galilei- invariant equations (1) can also 
be found in [3], [5], [19] -[23]. Galilean analogues of the Bargman-Wigner equations 
are presented in a recent paper [24]. However, these equations became incompatible 
whenever a minimal interaction with an external e.m. field had been introduced [24]. 

2 The Galilei algebra and Galilei— invariant wave 
equations 

2.1 Basic definitions 

In this section we shall develop a Galilean version of the Bhabha and Gel'fand-Yaglom 
approach [9] and present a complete list of the corresponding Galilei-invariant wave 
equations. Let us note that some of these equations are already known (see for 
instance, [3], [21] and [23]). 

Equation (1) is said to be invariant w.r.t. the Galilei transformations 

x -> x' = Rx + vt + b, t^t' = t + a (3) 

where a, b, v are real parameters and R is a rotation matrix, if function \1/ in (1) 
cotransforms as 

*(x, t) -> *'(x', if) = e i/(x ' * } T^(x, t), (4) 

i.e., according to a particular representation of the Galilei group. Here T is a matrix 
depending on transformation parameters only, /(x, t) = m (v ■ x + tv 2 /2 + c), c is 
an arbitrary constant and \l/'(x',t') satisfies the same equation in prime variables as 
\I/(x, t) in the initial ones. 

The Lie algebra corresponding to representation (4) has the following generators 

Po=Po = if, P a = Pa = -igjr, M = ml, 

Ja = SabcXbPc + S a (5) 
G a = tp a - mx a + T} a , 
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where S a and r] a are matrices which satisfy the following commutation relations: 

[Va,S b ] = ie abc r] c , [r} a ,r} b ] = 

that is, they form a basis of the homogeneous Galilei algebra hg(l, 3). 

Equation (1) is invariant with respect to the Galilei transformations (3), (4), 
if their generators (5) transform solutions of (1) into solutions. This requirement 
together with the existence of the Galilei-invariant Lagrangian for (1) yields the fol- 
lowing conditions on matrices (5^ (// = 0, 1, 2, 3,) and fa [22]: 

Tjlfa - far] a = -ifa, 

VlPb - PbVa = -iSabPo, (7) 

vlPo - PoVa = 0, a, b = 1, 2, 3. 

Moreover, fa and fa must be scalars w.r.t. rotations, i.e., they have to commute 
with S a . 

Thus the problem of classification of the Galilei-invariant equations (1) is equiva- 
lent to find matrices S a , rj a , fa, /3 a and fa satisfying relations (6) and (7). Unfortu- 
nately, a subproblem of this problem, i.e., a complete classification of non-equivalent 
finite-dimensional representations of algebra (6) appears to be in general unsolvable 
(that is a 'wild' algebraic problem). However, for two important particular cases, 
i.e., for purely spinor and vector-scalar representations, the problem of finding all 
finite-dimensional indecomposable representations of the algebra hg(l,3) is solvable 
and was completely solved in [5]. 



2.2 Spinor fields and the corresponding wave equations 

Let s be the highest value of spin which appears when representation of algebra 
hg(l,3) is reduced to its subalgebra so(3). Then the corresponding representation 
space of hg(l, 3) is said to be the space of fields of spin s. 

As mentioned in [5] there exist only two non-equivalent indecomposable repre- 
sentations of the algebra hg(l,3) defined on fields of spin 1/2. One of them, -Di(|), 
when restricted to the subalgebra so(3) remains irreducible while the other one, D 2 {\), 
decomposes to two irreducible representations .D(l/2) of so(3). The corresponding 
matrices S a and i] a can be written in the following form: 

Sa = \o a , Va = fOT D 1 (^j (8) 

and 

s ° = l(°o °)'"° = K° o) f ° rB2 G>)- 

Here a a are the Pauli matrices and is a 2 x 2 zero matrix. 
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Realization (8) with conditions (7) yields equation (1) trivial, i.e., with zero (3- 
matrices. 

The elements of the carrier space of representation (9) will be called the Galilean 
bi-spinors. It can be found in [5] how the Galilean bi-spinors transform w.r.t. the 
finite transformations from the Galilei group. 

Solutions of relations (7) with S a ,r) a given by formulae (9) can be written as 

*=(S o )•*.= (» o* )•*=(£ If) -= 1 ' 2 ' 3 ' < 10 > 

where I and are the 2x2 unit and zero matrices respectively, and u and k are 
constant multipliers. 

Notice that parameter k can be chosen zero since the transformation ^ — > e 1Kmt ^ 
leaves equation (1) invariant. Parameter uo is inessential too since it can be annulled 
by the transformation f3 m — > U^f3 m U, where m = 0, 1, 2, 3, 4 and 



U = 



I -icul 
I 



Let us note that if we consider a more general case in which matrices S a and r\ a 
are represented by a direct sum of an arbitrary finite number of matrices (9) and solve 
the related equations (7), then we obtain matrices (3 m which can be reduced to direct 
sums of matrices (10) and zero matrices. In other words, equation (1) with matrices 
(10) is the only non-decoupled system of the first order equations for spin 1/2 field 
invariant under the Galilei group. 

Equation (1) with matrices (10) and u = k = coincides with the Levy-Leblond 
equation in [3]. 

Let us remark that matrices % = r]P n \ K=OJ=0 , n = 0, 1, 2, 3, 4 with 

n-(°l) (ID 



x I o / 

satisfy the following relations 

7n7m + 7m7n = 2{jim, (12) 

where is a symmetric tensor whose non-zero components are 

g 4 = 940 = -hi = -922 = ~933 = 1- (13) 

In the Galilei-invariant approach tensor (13) plays the same role as the metric 
tensor (50) for the Minkovski space in the relativistic theory. 

The matrices % will be used many times later on. Therefore, for convenience, we 
present them explicitly, namely: 

7o = ( ° I ) , % = ( ° ~l a ) , 74 = ( o 2 q J, a = 1,2,3. (14) 
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2.3 Scalar and vector fields and the corresponding wave equa- 



tions 

2.3.1 The indecomposable representations for scalar and vector fields 

A complete description of indecomposable representations of the algebra hg(l,3) in 
the spaces of vector and scalar fields is given in [5]. The corresponding matrices S a 
and i] a have the following forms: 

s _( Inxn ® s a j ( A nxn ® S a B nxm ® k a 

\ " Omxm / \ Cmxn ® k a mxm 

where I nxn and mxm are unit and zero matrices of dimension n x n and m x m, 
respectively, A nxn , B nxm and C mxn are matrices of indicated dimensions whose forms 
will be specified later on, s a are matrices of spin one with elements (s a )bc = itabc and 
k a are 1x3 matrices of the form 

£4 = (i,0,0), fc 2 = (0, i, 0) , fc 3 = (0,0,i). (16) 

Matrices (15) fulfill relations (6), iff matrices A nxn , B nxm and C mxn satisfy the 
following relations (we have omitted the related subindices): 

AB = 0, CA = 0, A 2 + BC = 0. (17) 

This system of matrix equations appears to be completely solvable, i.e. it is possible 
to find all non-equivalent indecomposable matrices A, B and C which satisfy relations 
(17). Any set of such matrices generates a representation of the algebra hg(l, 3) whose 
basis elements are of the forms specified in (15). 

According to [5] indecomposable representations D(n,m,X) of hg(l,3) for scalar 
and vector fields are labelled by integers n, k and A. They specify dimensions of 
submatrices in (15) and the rank of matrix B, respectively. As shown in [5], there 
exist ten non-equivalent indecomposable representations D(n,k,X) of hg(l,3) which 
correspond to matrices A nxn , B nx ^ and Ckxn given in the Table 1. 

In addition to the scalar representation whose generators are written in (15) and 
in Item 1 of Table 1, there exist nine vector representations corresponding to matrices 
enumerated in Table 1, items 2-10. The corresponding basis elements are matrices of 
dimension (3n + k) x (3n + k) whose explicit forms are given in (15) and in Table 1. 
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Table 1. Solutions of equations (17) 



No 



I. 
2. 
3. 
4. 



(n, k, A) 



(0,1,0) 
(1,0,0) 
(1,1,0) 



Matrices^, B,C 



A, B and C do not exist since n = 

B and C do not exist since k = 

£ = 0, C= 1 
B = l, C = 

' 
1 




The finite Galilei transformations of vector fields (which can be obtained by in- 
tegrating the Lie equations for generators (15)) and examples of such fields can be 
found in paper [5]. 



2.3.2 General wave equations for vector and scalar fields 

Let us consider equation (1) and describe all admissible matrices fa compatible with 
the invariance conditions (7). We shall restrict ourselves to matrices r) a , S a belonging 
to the representations described in Subsection 2.3.1 (see Table 1) or to direct sums 
of these representations. Then the general form of matrices S a and r\ a is again given 
by equations (15) where, however, matrices A, B and C can be reducible: 



Mi 



A = 



\ 



( B 1 



B = 



\ 



B, 



c = 



Co 



\ 



■ / 

(18) 



and are of dimensions N x N, M x N and N x M respectively with N and M being 
arbitrary integers. The unit and zero matrices in the associated spin operator S 
defined by equation (15) are (N x N)— and (M x M)— dimensional, respectively. 
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The sets of matrices (Ax, Bi,C\), (A 2 , B 2 , C 2 ), ■■■ are supposed to be indecom- 
posable sets presented in Table 1. Any of them is labelled by a multi-index qi = 
(rij, fcj, Aj), 2 = 1, 2, • • • . 

Matrices /3 4 and /3 must commute with S and therefore have the following block 
diagonal form: 

RnxN OnxM \ a ( FnxN OnxM 



P _ j Jt -JVXiV "INX1V1 j p I J-iVXiV "INX1V1 j (19) 

Let us denote by \q, s, v > a vector belonging to a carrier space of the represen- 
tation D q of the algebra hg (1,3), where q = (n, k, A) is a multi-index which labels a 
particular indecomposable representation as indicated in Table 1, s is a spin quan- 
tum number which is equal to 0,1 and index v specifies degenerate subspaces with 
the same fixed s. Then taking into account that matrix fi± commutes with S a its 
elements can be written as 

< 9, s, v\f3 4 \q', s', v' >= 5 s i5 s >iR vu >(q, q') + 5 s0 5 sl0 E vu ,(q, q'). (20) 

In order to find matrices R(q, q') and E(q, q') (whose elements are denoted by R vv i 
and E vv i respectively) expression (20) has to be substituted into (7) and matrices r\ 
in form (15) together with relations (17) used. As a result we obtain the following 
condition 

(A ] fR + R(A'f = A ] RA' - C ] EC', (21) 

where A, C (and A', C) are submatrices used in (17), which correspond to repre- 
sentation D q (and D q i). 

Formulae (20) and (21) express all necessary and sufficient conditions for matrix 
f3& imposed by the Galilei invariance conditions (7). Suppose matrix (3^ in (20) which 
satisfies (21) be known, then the remaining matrices f3 a (a = 1,2,3) and /3 can be 
found by a direct use of the first and second relations in (7). By this way we obtain 

< q, s, \\PoW, s', A' >= 5 sl S sll F XX '(q,q') + <WW? A a'(<7, q'), 

< q, s, \\(3 a \q', s', A' >= i(5 sl 8 sfl H xy (q, q')s a + 5 sl 5 sl0 M xx ,(q, q')k\ (22) 

-5 sQ 8 sll Ml xi (q,q')k a ), 

where s a are matrices of spin one, k a are matrices (16) and F(q, q'), G(q, q'), H(q, q') 
and M(q, q') are matrices defined by the following relations 

H = A^R- RA', M = C*E - RB' , F = 0EC + A^RA', 

G = IB^RB' - B^C^E - EC'B'. ^ > 

Thus, in order to derive a Galilei- invariant equation (1) for vector fields, it is 
sufficient to choose a realization of the algebra hg(l, 3) from Table 1 or a direct sum 
of such realizations and find the associated matrix /3 4 (20) whose block matrices R and 
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E satisfy relations (21). Then the corresponding matrices fl and (3 a are determined 
via relations (22) and (23). 

All non-trivial solutions of matrices R and E are specified in the Appendix. Thus 
formulae (20)- (23) and Tables 2-4 in the Appendix determined all possible matrices 
(3± and f3^ which define the Galilei-invariant equations (1) for fields of spin 1. 

Notice that all these equations admit a Lagrangian formulation with Lagrangians 
in the following standard form 

L = -^ (^jf + 13 A m) ^ + h.c. (24) 
2.3.3 Consistency conditions 

In the previous section we have found all matrices fi m for which equation (1) is invariant 
with respect to vector and scalar representations of the homogeneous Galilei group. 
However, the Galilei invariance itself guarantees neither the consistency nor the right 
number of independent components. Moreover spin content of each of the obtained 
equations as well as their possibilities to describe fundamental quantum-mechanical 
systems have not yet been discussed . 

In this subsection we present further constraints on matrices f3 m which should be 
imposed in order to obtain a consistent equation for a system with a fixed spin. 

In this context, a non-relativistic quantum system is said to be fundamental if the 
space of its states forms a carrier space of an irreducible representation of the Galilei 
group (7(1,3). We shall call such systems as "non-relativistic particles" or simply 
"particles". 

A non-relativistic quantum system is said to be composed provided the space of 
its states forms a carrier space of some reducible representation of G(l, 3). 

First we consider equations (5) for fundamental systems. Notice that for the group 
(7(1,3) there exist the following three invariant operators: 

d — M, C 2 = 2MP - P 2 , and C 3 = (MJ - P x G) 2 . (25) 

Here, M and Pq are scalars and P, J and G are three- vectors whose components are 
specified by equations (5). Eigenvalues of the operators C±, C 2 and C3 are associated 
with mass, internal energy and with square of mass operator multiplied by eigenvalues 
of total spin, respectively. 

Galilei-invariant equation (1) is said to be consistent and describes a particle with 
mass to, internal energy e and spin s if it has non-trivial solutions \1/ which form a 
(non-degenerate) carrier space of some representation of the Galilei group on which 
the following conditions are true: 

Ci* = to*, C 2 ^ = and = m 2 s(s + 1)*. (26) 

Moreover, for fundamental particle the spin value s is fixed which corresponds to an 
irreducible representation, while for composed system the corresponding spin operator 
has two values: s = and s = 1 since the associated representation is reducible. 
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In accordance with the above definition of consistency, the number of independent 
components of function \1>, satisfying a consistent equation (1), must be equal to the 
number of spin degrees of freedom, namely to 2s + 1 for a fundamental particle and 
to 4 for a composed one. 

Now we shall show that relations (26) generate extra conditions for /^-matrices so 
that equation (1) with such (3 m guarantees the validity of equations (26). 

Function ^ satisfying equation (1) cotransforms according to a particular repre- 
sentation of the Galilei group whose infinitesimal generators are specified in (5) Using 
(5) we find the following forms of invariant operators (26): 



Ci = Im, C 2 = 2mp - p 2 , 

C 3 = m 2 S 2 + m(S x rf) ■ p - m(r] x S) ■ p + p 2 ry 2 - (p • r}) 2 



We see that C3 is a rather complicated second-order differential operator with 
matrix coefficients. In order to diagonalize this operator, we apply the similarity 
transformation 

_, = WV, C a ^C' a = WC a W-\ a = 1,2, 3, (28) 

with 

W = exp (j-rj ■ p^ . (29) 

Let us remark that since (77 • p) 3 = for representations D(3, 1, 1) and D(l, 2, 1) 
and (77 ■ p) 2 = for the remaining representations described in Subsection 2.3.1, W 
is the second- or the first-order differential operator in x. 

Using conditions (6) we find that 

C[ = d, C' 2 = C 2 and C'i = S 2 . (30) 

Function \P' has to satisfy conditions (26) with the transformed invariant operators 
(30), from which follow that 

(2m Po - p 2 )^' = e*' (31) 

and 

SV = s(s + (32) 

In order to see when conditions (31) and (32) are true we transform equation (1) 
by means of W. We obtain 

(p C 2 + (3 A 2m 2 ) = (33) 

since, in accordance with (7), we have: 

2m(W-y(^p II + fam)W- 1 = [3 (2mp - p 2 ) + (3 4 2m 2 . 
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Equation (33) in order to be compatible with (31) implies that the matrix (3 e + 
/9 4 2m 2 (where e is an eigenvalue of the Casimir operator C 2 ) should be non- regular 
for some particular value of e. Moreover, solutions of equation (33) must also satisfy 
condition (32) and form a carrier space of irreducible representation D{s) of the 
rotation group; consequently, equation (33) must have (2s + 1) independent solutions. 
Let us find now restrictions on matrices (3q and $4 generated by these solutions. 

As shown in Section 2.3.1, both matrices (3q and $4 have the block diagonal form 
given by equation (19). Thus equation (33) is decoupled to two subsystems 

(RC 2 + 2m 2 F) i Pl = 0, (34) 

and 

(EC 2 + 2m 2 G) v? 2 = 0, (35) 

where the functions (fii and ip 2 denote eigenfunction columns with 3n and k compo- 
nents respectively so that \P' = column (tpi, tp 2 ) . 

Equations (31) and (32) are consequences of (34) and (35) provided determinants 
of matrices R, F, E and G satisfy the following conditions: 

det (RC 2 + 2m 2 F) = vC 2 + det [EC 2 + 2m 2 G) = v'C 2 + //, (36) 

where //, v and //, v' are numbers satisfying one of the following relations: 

i/^0, fi= -ev, v' = 0, for s = 1, (37) 

u = 0, fj, 0, v' ^ 0, n' = -ev' for s = 0. (38) 

Thus to find Galilei- invariant equations (1) for particle with a fixed spin we can 
start with equations described in the previous section and impose conditions (36) on 
block components of matrices (3q and $4. 

In analogous way we find the consistency conditions for equations (1) describing 
composite systems. In this case relations (31) and (32) should be valid too. However, 
s can take two values s — 1 and s = 0. We suppose that these spin states are non- 
degenerate, so that function \J/' should have four independent components - three 
of them corresponding to spin one and one to spin zero state. Then conditions for 
determinants of matrices FC 2 + 2m 2 R and GC 2 + 2m 2 E again have form (36) where, 
however, 

vv' ^ 0, pi — —ev, /i' = —ev'. (39) 

In the next sections we consider examples of various Galilean equations (1) in 
more detail. 
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3 Special classes of the Galilean wave equations 
for free fields 



3.1 Equations invariant with respect to the indecomposable 
representations 

In Sections 2.3.2 and 2.3.3 we have given a general description of all possible Galilei 
wave equations (1) for vector and scalar fields. Now we shall present an analysis of 
these equations in details. In this and next sections we shall consider equations for 
systems with both one fixed spin state and two spin values states. 

First we restrict ourselves to the indecomposable representations of the algebra 
hg(l, 3) specified in equation (15) and Table 1, and find the associated matrices Ao 
and P a which appear in the Galilei-invariant equations (1). Taking into account that 
A' = A, B' — B and C = C in (21), where A, B and C are matrices given in 
Table 1, and using the results present in Tables 2-4 in the Appendix, we easily find 
the associated block matrices R, E and consequently all block matrices (23) and 
matrices (20) and (22). In order to simplify matrices Am (m = 0, 1, • • • ,4) we use the 
transformations 



Am -> V%V, 



(40) 



where V are invertible matrices commuting with the Galilei boost generators r] a . 

It turns out that non-trivial solutions for /3 4 (and consequently also for Ao and 
P a ) exist only for representations D(l, 1, 0), D(2, 1, 0), D(2, 2, 1) and D(3, 1, 1). They 
have the form: 



For representation D(1,1,0): 



Pa 



13x3 
0lx3 



3 xl 





Ao 



03x3 
0lx3 



For representation D(2, 1,0): 



Pa 



Pa 



3 x3 


3 x3 


3 xl 


3 x3 


13x3 


3 xl 


0lx3 


0lx3 




3 x3 


Sa 




~S a 


3 x3 


3 xl 




0lx3 






Po 



03x1 
2 



Pa 



03x3 

kn 








(41) 



(42) 



For representation .0(2,2,1): 

/ 03x3 O3X3 O3XI O3XI \ 
03x3 13x3 O3XI O3XI 
0lx3 0ix3 

\ o lx3 o lx3 1 / 



Pa = 



, Po 



( 


213x3 


3 x3 


3 xl 


3 xl 


\ 




3 x3 


3 x3 


3 xl 


Osxl 






0lx3 


0lx3 


2 







V 


0lx3 


0lx3 








/ 
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/ 3X3 


Sa 


3 xl 


kl \ 


-So 


3 x3 


-kl 


3 xl 


0lx3 


k a 










0lx3 








Pa 



And finally for representation D(3, 1, 1) : 



(43) 



Pa 



Pa 



03x3 
^3x3 
^3x3 
0lx3 

3 x3 
3 x3 
03x3 

k a 



^3x3 
x 3x3 
3 x3 
0ix3 



03x1 
3 xl 
3 xl 
— V 



\ 



Po 



J 



( o 3x3 


13x3 


03x3 


3 xl 


\ 


13x3 


o 3x3 


3 x3 


3 xl 




3 x3 


3x3 


3 x3 


3 xl 




V 0lx3 


0lx3 


0lx3 





/ 



Sa 
3 x3 
03x3 
0lx3 



3 xl \ 
3x i 

/ 



(44) 



/ 3x3 

o 3x3 

^3x3 

V o lx3 

/ o 3x3 
o 3x3 

\ 0lx 3 

Here, Ifc X fc and 0k Xr are unit and zero matrices of dimensions k x k and k x r,respectively 
and v is an arbitrary non-vanishing parameter. 

Thus there are four equations (1) for spinor and vector fields which are invariant 
with respect to the above mentioned indecomposable representations of the homo- 
geneous Galilei group. Their associated matrices P^ and P4 are given by formulae 
ill). 

The equations (1) with /3-matrices specified in (41), (42), and (44) are equivalent 
to those discussed in papers [20] and [21], [22] respectively. However, equation (1) 
with /3-matrices of the form (43) is to the best of our knowledge new. Notice that 
the related submatrices F, R and G, E satisfy relations (36) and (39), so the associate 
equation describes a Galilean quantum-mechanical system whose spin can take two 
values: s = 1 and s = 0. 

Considering matrices (44) we conclude that the corresponding matrices 



Pn = VPv, V- = 0, 1, 2, 3, and /3 4 = t]Pa 
where 77 is an invertible matrix 



V = 



satisfying the following relations: 

PuPvPv + PaPvPp, = ^{g^uPa + QauPfi), 



IOxIOi 



(45) 



^ 03x3 


03x3 


^3x3 


O3XI 


\ 


03x3 


^3x3 


03x3 


O3XI 




^3x3 


03x3 


03x3 


03x1 




V o lx3 


0lx3 


0lx3 


-1 





(46) 



(47) 



where is a "galilean metric tensor" given in (13). 

Up to definition of metric tensor relations (47) coincide with the defining 
relations for the Duffin-Kemmer-Petiau (DKP) algebra [25]. Following [26] we say 
that such relations define a Galilean DKP algebra. 
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Equation (1) with matrices fa, fa which satisfy the Galilean DKP algebra is called 
the Galilean Duffin-Kemmer equation. This equation was considered for the first time 
apparently in [21]. 

There exist also a number of wave equations (1) invariant with respect to decom- 
posable representations of hg(l, 3). We shall discuss some of them in the subsections 
which follow. 

3.2 The Galilean Proca equations 

In addition to equations discussed in the previous section (whose solutions belong to 
carrier spaces of indecomposable representations of the homogeneous Galilei algebra), 
there exist also a number of wave equations (1) invariant with respect to decomposable 
representations of hg(l, 3). Complete description of such equations is given in Section 
2.3.2 and the Appendix A. In this section we shall discuss some of them as examples. 

Let us consider first wave equations whose solutions are vectors from a carrier 
space of the direct sum of representations D(3, 1, 1)©D(2, 1, 1) of the algebra hg(l, 3). 
The associated matrices fa, fa and fa can be found using relations (20), (22), (23) 
and Tables 2 and 3 where submatrices R(q, q), R(q, q'), R(q', q') and E(q, q), E(q, q'), 
E(q', q') are specified for q = (3, 1, 1) and q' = (1, 2, 1). We can simplify these matrices 
by using the equivalence transformations (40). 

In this way we realize that there are two consistent equations for the representa- 
tion -D(3, 1, 1) © -D(2, 1, 1): one describing a particle with spin s = 1 and the other 
describing a composed system with spins s = 1 and s = 0. Here we present a covariant 
formulation of these equations. 

The equation for a particle with spin s = 1 can be written in the following form 



Here, <5 nk is the Kronecker symbol, \& nk (n,k = 0,1,2,3,4) is an antisymmetrical 
tensor and \l/ n is a five-vector which transform in accordance with representations 
D(3, 1, 1) and D(l, 2, 1), respectively. 

Equation (48) includes an arbitrary real parameter A ^ whose value cannot be 
fixed using only conditions of the Galilei invariance. 

The equation for a system with two spin states involves again tensor and five- 
vector variables which satisfy the following system: 



The Galilei-invariant equations (48), (49) are quite similar to the relativistic Proca 
equations (2), but with the following differences: 

• all indices m,n in (48) take values 0,...,4 while in (2) we have /i, v — 0, 1, 2, 3; 



= A<5 n4 m^ 4 . 



(48) 



j9 k \J/ n _ jr, n \J/ k — m \fy 



(49) 
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the relativistic four-gradient p v is replaced by the Galilean five-vector p = 

■a_ a _ 



(p°, p 1 , p 2 , p 3 , p 4 ), where p° = i4r, p a = ijf-, a = 1, 2, 3, and p 4 = m; 



• covariant indices m, n are raised and lowered by using the Galilean metric tensor 
(13) instead of the relativistic metric tensor 

g llv = diag(l, -1,-1,-1); (50) 



• finally, equation (49) describes a system with two spin states, while both the 
relativistic Proca equation and the Galilei-invariant equation (48) describe a 
particle with fixed spin s — 1. We shall call (48) the Galilean first order Proca 
equation. 

Let us stress here that the considered Galilean Proca equations are not non- 
relativistic approximations of the relativistic Proca equation (2) (see appendix C 
where a contraction of (2) to its non-relativistic approximation is shown). 

Substituting the expression for \l/ nk from the first into the second equation in (48), 
we obtain the Galilean second order Proca equation: 

p n p n ty m - p m p n ty n + X5 m °m 2 * 4 = (51) 

where m,n = 0,1,2,3,4. Equation (51) admits, like (48), a Lagrangian formulation 
and describes a particle with spin 1. The corresponding Lagrangian has the form: 

-(pA)y^ n -Am 2 $^ 4 , 1 ' 

where the asterisk * denotes complex conjugation. 



3.3 The Galilean Rarita-Schwinger equation 

Till now we have used our knowledge of the indecomposable representations of the 
algebra hg(l, 3) for spinor, scalar and vector fields to construct wave equations for 
fields of spin s < 1. In this section we derive Galilean invariant equations for the 
field transforming as a direct product of spin 1/2 and spin 1 fields. The relativistic 
analogue of such system is the famous Rarita-Schwinger equation. 

The relativistic Rarita-Schwinger equation for a particle with spin s — | is 
constructed by using a vector-spinor wave function where p, = 0,1,2,3 and 
a = 1,2,3,4 are vector and spinor indices respectively. Moreover, ty£ is supposed 
to satisfy the equation 

(YPu -m)W- ^p v W - p^ v + Y { lv p v + m) 7a* CT = 0, (53) 

where 7 M are the Dirac matrices acting on the spinor index a of which we have 
omitted. 
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Reducing the left-hand side of equation (53) by p M and 7 M we obtain the following 
expressions: 

7„*" = 0, and = 0, (54) 

which reduce the number of independent components of to 8 as required for a 
wave function of a relativistic particle with spin 3/2. 

Using our knowledge of invariants for the Galilean vector fields from [5] and [6] we 
can easily find a Galilean analogue of equation (53). Like in the case of the Galilean 
Proca equation we begin with a five- vector \I/ m , m = 0, 1, 2, 3, 4 which has, in addition, 
a bi-spinorial index which we do not write explicitly. Thus our Galilei-invariant 
equation can be written in the following form: 

7n p n £ m - 7 > n ^ n - p m 7n ^ n + 7 m 7 n p n 7 r ^ r + A5 m( W 4 = 0. (55) 

Here, % are the Galilean 7-matrices (14), p m is a Galilean "five- momentum" defined 
in the previous subsection (see the second Item there), A is an arbitrary non- vanishing 
parameter and raising and lowering of indices m and n is made by using the Galilean 
metric tensor (13). 

Like equation (51) the equation (55) admits a Lagrangian formulation. The cor- 
responding Lagrangian can be written as 

L = |(^ m 7nP n ^ m - ^ m 7>n^ n ~ ^^ln^ + ^mTlnP^la^ + Am\^ 4 ) + h.C, 

(56) 

where \fr m = and i] is the hermitizing matrix (11). 

We prove in Appendix B that equation (55) indeed describe a particle with spin 
3/2. 



4 Equations for charged particles interacting with 
an external gauge field 

4.1 Minimal interaction with an external field 

We have described Galilei-invariant equations (1) for free particles with spins 0, 1/2,1 
and 3/2. These equations have admitted Lagrangian formulation (24), so that to 
generalize them to the case of particles interacting with an external field means, as 
usually, to apply the minimal interaction principle, i.e., to make the following change 
in the Lagrangian 

p" -> tt" = - eA», (57) 

where are components of a vector-potential of the external field, and e is a particle 
charge. 
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Thus we obtain the Lagrangian 

L = ltft( x> t) + Am) *(x, + h.c. (58) 

It is important to note that change (57) is compatible with the Galilei invari- 
ance provided the components (A , A 1 , A 2 , A 3 ) of the vector-potential transform as a 
Galilean four- vector, i.e., as 

A°^A° + v-A, A^A. (59) 

Such potential corresponds to a "magnetic" limit of the Maxwell equations, see [19], 
[27]. For other possible Galilei potentials of the electromagnetic field see [7]. 

The desired equations for a charged particle interacting with an external field are 
the Euler equations derived from the Lagrangian (58). 

However let us note that it is also possible to introduce interaction with an external 
field via other means, e.g., via an anomalous (Pauli) term. 



4.2 The Galilean Bhabha equations with minimal and anoma- 
lous interactions 

Taking Lagrangian (58) we can derive the following equation for a charged particle 
interacting with an external field 

(/3X + &7r 4 )tt(x,i) = 0. (60) 

Let us consider now equation (60) with general matrices (3^ and (3±. If we restrict 
ourselves to a vector-potential of magnetic type, i.e., to A = (A , A, 0), then 

n° =p° - eA°, n a =p a - eA a , tt 4 = m. (61) 

Like free particle equations (1), the equation (60) is Galilei-invariant provided 
matrices #i satisfy conditions (7). Moreover the vector-potential of an external 
field has to transform according to (59). 

Following Pauli [29] we generalize our equation (60) by adding to it an interaction 
terms linear in an electromagnetic field strength. Then we get the equation: 

(/Vr" + Am + F) = 0, (62) 

where 

F = — (A-H + G -E). 

m 

Here, A and G are matrices determined by requirement of the Galilei invariance, i.e., 
by demanding that A • H and G • E have to be Galilean scalars. 
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In paper [5] we have found the most general form of the Pauli interaction which 
can be introduced into the Levy-Leblond equation for a particle of spin 1/2 . Finding 
the general Pauli interaction for other Galilean particles is a special problem for any 
equation previously considered. Here we restrict ourselves to a systematic analysis of 
the Pauli terms which is valid for any Galilean Bhabha equation. 

First we shall prove the following statement. 

Lemma. Let S a ,r) a be matrices which realize a representation of the algebra 
hg(l,3), A be a matrix satisfying the conditions 

S a A = ^S a , vi^ = ^Va (63) 

and 

_dAo_dAa „ _ dA^ 
a dx- dt ' a £abc dx c 

be vectors of the electric and magnetic field strength, respectively. Then matrices 

Fi — A(s • H — 77 • E) and F 2 = Arj ■ H (64) 

are invariant with respect to the Galilei transformations provided the vector-potential 
A is transformed in accordance with the Galilean law (59). 

Proof. First we note that matrices (64) are scalars with respect to rotations. 
Then, starting with transformation laws (3) and (59) we easily find that under a 
Galilei boost the vectors E and H co-transform as 



E^E-vxH, H^H. (65) 

On the other hand the transformation laws for matrices AS and At] can be found 
using the exponential mapping of boost generators G given in equation (5): 

S -> exp(iGt • v)AS exp(-iG • v) 
= A exp(ir) • v )S exp(— irj • v ) = s + v x 77, , , 

r] — > exp(iG^ • v)Ar7exp(?G • v) 
= Aexp(i/7 • v )rf exp(-ir) • v ) = Ary . 

One easily verifies that transformations (65) and (66) leave matrices Fi and F 2 
invariant. Q.E.D. 

In accordance with the Lemma there are many possibilities how to generalize equa- 
tions (60) to the case with anomalous interaction. Indeed, for any Galilean Bhabha 
equation there are matrices S a , i] a and A for which conditions (63) are satisfied. For 
example, we can choose A = /3q. In addition, for many cases there exist a hermitiz- 
ing matrix rj = A satisfying (63), see, e.g., equations (11) and (46). For particular 
representations of the algebra hg(l, 3) there are also other solutions of equations (63). 
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Thus the Pauli term for a Galilean Bhabha equation can be chosen in the form 
(64) or, more generally, as a linear combination of both, F\ and F2. As a result we 
obtain the following equation, 

Q* = (/3„7r" + fam + X 1 ^ oV • H + A 2 ^ (S • H - 77 • E)) * = 0, (67) 

where Ai and A2 are dimensionless coupling constants. 

Since equation (67) is reduced to equation (60) for Ai = A2 = 0, equation (67) 
describes anomalous as well as minimal interaction. 

In order to receive the physical content of this equation it is convenient to apply 
the transformation 

* -> = W'H, Q^Q' = W^QW, (68) 

where 

W = exp (-v^-^) (69) 

and 77 is a vector whose components are the Galilei boost generators (9). For the 
case e = (or = 0) the operator W reduces to operator U given in equation (29), 
which was used for our analysis of free particle equations. 

Using relations (7) and supposing that the nilpotence index N of matrix r? ■ 7r 
satisfies iV < 4 we obtain the following equation: 



An, .H + A 2 (s.H-„ -F + ^Q ab ^)]}r = 



(70) 



which is equivalent to (67). Here E = — VA° — ^ and H = V x A are vectors of the 
corresponding electric and magnetic field strength respectively, 

F = E+2^(tt xH-Hxtt) (71) 

and 

Qab = ^{VaSb + VbS a + S b i] a + i] a S b ), Q ab = + e bcd (3 c i] d i] a . 

Equation (70) includes the Schrodinger terms (7r° — ^)^' and additional terms 
which are linear in vectors of the external field strengthes and their derivatives. 

Notice that if the nilpotence index N of matrices i] a satisfies the relation N < 4, 
which is fulfilled for all representations of algebra hg(l,3) considered in the present 
paper, the transformed equation (70) is completely equivalent to initial equation (67). 
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4.3 The Galilean equation for spinor particle interacting with 
an external field 

Let us consider equation (70) for two particular realizations of /3-matrices in more 
detail. First notice, that our conclusions from equations (60)-(70) are true in general 
and in particular for the Levy-Leblond equation, i.e., when (3^,04 are 4x4 matrices 
determined by relations (10) with u — k — 0. Then (3 rj a = 0, Q a b = Q a b = 0, f3xrj = 
—2j3 S, and equation (70) is reduced to the following form: 

{A^-S + ^-Hj+Am + ^AfAn, H + A 2 (S H-r7 -F)]}*' = 0. 

(72) 

For Ai = A 2 = (i.e., when only the minimal interaction is present) equation (72) 
is reduced to the following system: 

*>-£ + £r- H )'" = (73) 

and 

rmp 2 = 0, or ip 2 = 0, (74) 

where ipi = (3q^' and <f 2 — (1 — Po)^' are two-component spinors. 

Thus introducing the minimal interaction (57) into the Levy-Leblond equation, 
we get the Pauli equation for physical components of the wave function; moreover, 
the coupling constant (gyromagnetic ratio) for the Pauli interaction • H, s = \a 
has the same value g = 2 as in the case of the Dirac equation [3]. 

Considering now the case with an anomalous interaction we conclude that the 
general form of matrix A satisfying relations (63) is 

A = vfio + fir], (75) 

where 77 is hermitizing matrix (11), v and /i are arbitrary parameters. Substituting 
(75) into (72) we obtain the following generalization of system (73): 



7T 



£ + ^-h-^'F-Sh 2 ) ^ = 0, 



V2 = -S"-H¥>i, 



(76) 



where g = 2 + //Ai + v\ 2 and A3 = fi\ 2 are arbitrary parameters. 

We see that the Levy-Leblond equation with minimal and anomalous interactions 
reduces to the Galilean Schrodinger-Pauli equation (76) which, however, includes two 
additional terms linear in strengths of an electric field and linear and quadratic in 
strengths of a magnetic field. We shall discuss them in detail in Section 4.6. 
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4.4 The Galilean Duffin-Kemmer equation for particle inter- 
acting with an external field 

Consider now the Galilean Duffin-Kemmer equation for spin one particles interacting 
with an external field. The corresponding /^-matrices in (67) have dimension 10 x 10 
and are given explicitly in (44). Thus 



( S 03x3 03x3 

03x3 03x3 03 X 3 

03x3 03x3 03x3 

\ 0ix3 0ix3 0ix3 



Q 



ab 



( —"iQab 
03x3 
03x3 

V o lx3 



03x3 
03x3 
03x3 
0lx3 



03x3 
03x3 
03x3 
0lx3 



3 xl \ 
3 xl 
03x1 

/ 

3 xi \ 
03xi 
03xi 




A)S 



^03x3 


S 


03x3 


03x1 


\ 


S 


03x3 


03x3 


03x1 




03x3 


03x3 


03x3 


03x1 




V 0ix3 


0lx3 


0lx3 





/ 



, Q 



ab 



( Qab 3X 3 

03x3 03x3 

03x3 03x3 

\ 0ix3 0ix3 



03x3 
03x3 
03x3 
0lx3 



(3 xr/ H = - 



/ 


3X3 


s H 


03x3 


2kt • H \ 




s H 


03x3 


03x3 


O3XI 




03x3 


03x3 


03x3 


03x1 


V 


2k H 


0lx3 


0lx3 


/ 



3 xi \ 
03xi 
03xi 
/ 



(77) 



where 



Qab S a Sfe 4" S b S a 



-5 ab and Q' ab = Q ab + -5 ab . 



(78) 



Let us consider the corresponding equation (70) and restrict ourselves to A = (3 . Rep- 
resenting iff' as a column vector (fa, fa, fa, V 9 )) where fa, fa, fa are three-component 
vector functions and ip is a one-component scalar function and using (44), (77), (78) 
we reduce (70) to the following Pauli-type equation for fa: 



Hfa, 



where 



H 



T-+l^ + ^o-||s.H + ^s.E- 
lT^(2-^ff + ^(H 2 -(s.H) 2 ) 



2um 2 



(79) 



S • (-7T X H - H X 7r) 



(80) 



and g — 1 + 2Ai + 2A 2 , q = 1 — A 2 . 

The remaining components of iff' can be expressed in terms of fa: 



? k -Ufa, fa = -\fa, fa = -Ufa - ± (tT - ^7T 2 + • H) fa. 
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Notice, that in comparison with (76) the equation (79) has an essentially new 
feature. Namely, excluding anomalous interaction, i.e., setting Ai = A 2 = in (80) 
it still includes the term — • E describing the coupling of spin with an electric 
field. We shall show in the next section that this effectively represents the spin-orbit 
coupling. The other terms of Hamiltonian (80) (which are placed in the second line of 
equation (80)) can be neglected starting with a reasonable assumption about possible 
values of the magnetic field strength. 

However, equation (76) for Ai = A 2 = reduces to the Schrodinger-Pauli equation 
(73) which has nothing to do with the spin-orbit coupling. 



4.5 The Galilean Proca equation for interacting particles 

Equations (48) and (49) can be generalized too by introducing minimal and anomalous 
interactions with an external field. As a result we obtain: 



7r k # nk = um^ n + i/i^F nk ^ k + A5 n4 m# 4 , 



(81) 



where parameters A and v satisfy the conditions \v = 0, u 2 + A 2 ^ 0. Formulae (81) 
generalize both equation (48) (which corresponds to v = 0, A ^ 0) and equation (49) 
(for which v ^ 0, A = 0). 

Multiplying equations (81) by both 7T n and 7r n 7r k , summing up by m and k and 
then expressing \l/ nk , \1/ via \l/ a and \1/ 4 , we obtain the following system: 

(To-£ + 1 S i s.H-|m)* + ^(l-^)kt.F^ = 0, 

(82) 

(^o + ^)* 4 -^(l-^)k-F* = 0, 

where we have denoted ^ = column^ 1 , \P 3 ). 

Let v — then solving the second of equations (82) for \1/ 4 and substituting the 
result into the first equation we obtain 

(to - £ + (1 + 2/i)^s • H - |m + (1 - /i) 2 ^(F 2 - (s • F) 2 ) ¥ = 0. (83) 

In accordance with (83) the Galilean Proca equation for a particle with spin 1 
interacting with an external field is reduced to the Schrodinger-Pauli equation with 
an extra term ^-(F 2 — (s • F) 2 ) which can be treated as a small correction. 

For v 7^ equations (82) describe a composed system with spins s = 1 and s = 0. 
Notice that there are two privileged values of the arbitrary parameter u, namely, 
v — — 1 and v — 1. 

If v — — 1 or v — 1 then equations (82) can be rewritten in the Schrodinger form 
(79), where ^i = column^ 1 , ^ 2 , ^ 4 ), and 

H = f + eA -^S-U- 3 —^K-F-^ (84) 
2m 2m 2m 2 v ; 
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or 

71-2 a Q£ „ TT 3 — q - ,-, ?7i ,„_ s 
# = + eA -— S-H -K-F + — , 85 

2m 2m 2m 2' v ; 

respectively. Here g = 1 + 2fi and S, K K are matrix vectors whose components are 

3x i \ J, ( 3x3 K\ T > ( °3xs kt 



s - = { oL o ' J ■ - { T o ) • - { -Z J < 86 > 

and k a are matrices (16). Matrices {S a ,K a } and {S^, lf } form bases of the algebra 
so(4) and so(l,3), respectively. 

Notice that Hamiltonian (84) is formally Hermitian w.r.t. the standard scalar 
product for the direct sum of four square integrable functions while (85) is Hermitian 
w.r.t. the indefinite metric (^i, fa) = f ^{MfacPx, where M is either a matrix whose 
elements are given in equation (50) or a parity operator. For other non-vanishing 
values of the parameter v in (81), i.e., for v ^ 0, ±1, the corresponding Hamiltonian 
H is non-Hermitian. 



4.6 The Galilei invariance and spin-orbit coupling 

Consider now the first of equations (76) for particular values of arbitrary parameters, 
namely for Ai + A 2 = — 1: 

M--^-H-' f -IH^ (87) 

First, let us remind that this equation is a direct consequence of the Galilei- 
invariant Levy-Leblond equation with an anomalous interaction, i.e., of equation (67) 
where /3 n are matrices (10) with k — uj — 0. Secondly, equation (87) by itself is trans- 
parently Galilei-invariant since the operator L in (87) is a Galilean scalar provided 
the value of an arbitrary parameter A3 is finite. We shall assume A3 to be small. 

In order to find out the physical content of equation (87) we transform it to a more 
transparent form using the operator U = exp(— ^-cr • n). Applying this operator to 
ipi and transforming L — > L' — ULU^ 1 we obtain the equation 

LV'i = (tto - m - eA - |J (<r • (tt x E - E x tt) - divE) + • • •) <p[ = 0, (88) 

where the dots denote small terms of the order o(A|) and o(e 2 ). 

All terms in big round brackets have an exact physical meaning. They include 
first the Schrodinger terms 7T — j- — eAo, then the term ~s-(7rxE-Ex7r) 
describing a spin-orbit coupling and, finally, a term ~ divE, i.e., a Darwin coupling. 

Similarly, starting with equation (79), setting A 2 = —1, Ai = |, supposing ^ to 
be a small parameter and making use of the transformation fa — > fa x —Ufa with 
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U = exp(— -^s ■ 77 ) we obtain the equation 

(*o - %m- £ -eA - £ (, • (77 x E-E x 77) + fdivE 

-^ff)+-->i = 0- 

Here, A = and the dots denote small terms of the order 0(^3) and o(e 2 ). 

Like equation (88), the equation (89) includes the terms which describe the spin- 
orbit and Darwin couplings. In addition, there is the term ~ Qaftff^ which describes 
a quadrupole interaction of a charged vector particle with an electric field. 

Let us stress that these terms are kept for the minimal coupling also, i.e., when 
we set Ai = A 2 = in (67) and (80). 

Analogously we can analyze equations (79) with the Hamiltonias (84) and (85). 
Setting there g — and making transformations H — > UHU" 1 — ill^—, where 
U = exp(— 3iK • 7z/2m) and U = exp(— 3iK • 77 /2m) for Hamiltonian (84) and (85), 
respectively, we obtain the approximate equation (89) with A = 9z//8, v = =pl. 

Thus we again come to the conclusion (see [21]) that the spin-orbit and Darwin 
couplings can be effectively described within the framework of a Galilei-invariant 
approach and thus they have not be necessarily interpreted as purely relativistic 
effects. 

Let us note that it is possible to choose parameters Ai and A2 in (80) in such a 
way that the anomalous interaction with an electric field will not be present. Namely, 
we can set A 2 = 1, Ai = —1/2, and obtain, instead of (89), the following equation : 

(^o-^m-g-eA) + ||s-H + --.)v>i = 0, (90) 

where g = 2. In other words, introducing a specific anomalous interaction into the 
Galilean Duffin-Kemmer equation we can reduce it to the Schrodinger-Pauli equation 
with the correct value of gyromagnetic ratio g. 



5 Discussion 

It is pretty well known that the correct definition of non-relativistic limit of relativistic 
theories is by no means a simple problem. In particular, as it was noticed once more 
in the recent paper [32], a straightforward non-relativistic expansion in terms oi v/c 
leads to loosing either Galilei-invar iance or important contributions such as spin- 
orbit coupling. Thus to make this limit correctly it is imperative to have a priory 
information on possible Galilean limits of a given theory. 

In the present paper we continue the study of the Galilei-invariant theories for 
vector and spinor fields, started in [5]- [7]. The peculiarity of our approach is that, 
as distinct to the other approaches (e.g., to [19]- [23], [28] ), it enables to find out a 
complete list of the Galilei-invariant equations for massive scalar and vector fields. 
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This possibility is due to our knowledge of all non-equivalent indecomposable repre- 
sentations of the Galilei algebra hg(l,3) that can be constructed on representation 
spaces of scalar and vector fields, i.e., the representations which were described for 
the first time in paper [5]. Notice that Galilean equations for massless vector and 
scalar fields are presented in our previous paper [7]. 

Using this complete list of representations we find all systems of the first order 
Galilei-invariant wave equations (1) for scalar and vector fields. The /5-matrices for 
these Galilei-invariant wave equations are given in the Appendix A. In fact we have 
described how to construct any wave equation of finite order invariant with respect 
to the Galilei group since such equation can be written as a first order differential 
equation in which various derivatives of fields are considered as new variables of 

Then Galilean analogues of some popular relativistic equations for vector particles 
and particles with spin 3/2 are described, in particular, the Galilean Proca and the 
Galilean Rarita-Schwinger equations. However these Galilean equations are not non- 
relativistic limits of the relativistic Proca or relativistic Rarita-Schwinger equations 
since, among other things, they have more components. Thanks to that it is possible 
to obtain equations which keep all the main features of their relativistic analogues. 
To the best of our knowledge this is done for the first time in the present paper. 

We pay a special attention to description of the Galilean particles interacting with 
an external electromagnetic field. We study both the cases with a minimal interaction 
as well as anomalous one. 

A quite general form of an anomalous interaction which satisfies the Galilei in- 
variance condition is written in equation (67). It contains two coupling constants, Ai 
and A 2 , whose values can be fixed via physical reasonings. If we fix the value of gy- 
romagnetic ratio predicted by (67) then the randomness in description of anomalous 
interaction is reduced to one arbitrary parameter. It can be fixed too if we restrict 
ourselves to a desired value of the spin-orbit coupling constant. 

Notice that the results presented in Sections 4.1 and 4.3 are valid for arbitrary 
equations (67) invariant with respect to the Galilei group whereas those in Sections 
4.2, 4.4 and 4.5 are true for special equations with anomalous interactions, i.e., for the 
Galilean Proca equation, the generalized Levy-Leblond equation and for a generalized 
Galilean Duffin-Kemmer equation. We shall show that the last equations describe 
consistently charged particles interacting with an electromagnetic field. In other 
words, they describe an important physical effect, namely, spin-orbit coupling which 
is, however, traditionally interpreted as a purely relativistic phenomenon. 

On the other hand let us note that there are some principal difficulties with the 
Galilean approach since the Galilei invariance requires that mass and energy are 
separately conserved, and that within the Galilean theories there is not concept of 
proper time which yields a phase effects that does not depend on the velocity of light 
and so does not disappear in a non-relativistic limit. Of course, there are obvious 
restrictions to phenomena which are characterized by velocities much smaller than 
that of light. Moreover, there are also problems in our approach with interpretation of 
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undesired terms ~ Qab^f and s • (jr x H — H x tz) which appear in Hamiltonian (79). 
Thanks to an appropriate choice of otherwise arbitrary parameters Ai and A 2 these 
terms are not present in effective Hamiltonians (89) as well as (90) which describe 
spin-orbit and the Pauli couplings respectively. However, if we would like to keep 
both these couplings, then the undesired terms may appear. 

Another problem is connected with a sign in front of the terms describing the 
spin-orbit and the Darwin couplings. Comparing (88) with the quasi-relativistic ap- 
proximation of the Dirac equation we conclude, that in order to obtain a correct 
signs it is necessary to suppose that A3 be pure imaginary. Moreover, for A3 = i the 
coupling constants for spin-orbit and the Darwin interactions in (88) coincide with 
the relativistic ones predicted by the Dirac equation. 

Notice that the exact equation (87) is much simpler then the approximate equation 
(88) and can be solved exactly for some particular external fields (for instance, the 
Coulomb ones). However, if A 3 is imaginary, the term — |^<r ■ F in equation (87) 
and the corresponding Hamiltonian H = —L + po = Aq are non-hermitian. If Aq 
and A are even and odd functions of x, respectively, the equation (87) appears to 
be invariant with respect to a product of the space inversion P and the Wigner time 
inversion T (compare with [33]), and thus serves as an example of a PT - symmetric 
quantum-mechanical system. 

Any Galilean theory, by definition, is only an approximation of the corresponding 
relativistic one. Thus the very existence of a physically consistent non-relativistic 
approximation can serve as a consistency criterium of a relativistic theory. Thus 
our study of the Galilean wave equations has contributed to the theory of relativistic 
equations, since effectively we have analyzed possible non-relativistic limits of theories 
for vector and scalar particles. On the other hand the results of the present paper can 
be applied to purely non-relativistic models satisfying the Galilei invariance criteria 
or being invariant w.r.t. various extensions of the Galilei group, e.g., w.r.t. Galilei 
supergroup [34]. 

A Submatrices R and E of matrices $4 

Here we present all non-trivial solutions of equations (21) which give rise to explicit 
forms of matrices ^4 given by equation (20). The associated matrices (5q and /3 a are 
given by equations (22) and (23). 

Solving equations (21), where A, C and A', C are matrices given in Table 1 which 
correspond to q — (n, m, A) and q' = (n f , m', A') respectively we obtain the associated 
matrices R = R(q, q'),E = E(q, q') (which determine matrix (3 4 via (20)) in the forms 
presented in Tables 2-4, where the Greek letters denote arbitrary real parameters. 
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Table 2. Submatrices R and E of matrices (3± 



m', n', X' 


to, n, A 


3,1,1 


2,2,1 


2,1,0 


3,1,1 


/ n v a \ 
R= \ v a A 
\ a X / 

E = a — 2(7, 

yuz/ = 0, X(a — a) = 


R =(" ° n ) 
^i/ a y 

r ^ i 

y cj — ct y 


y f a a; J 

E = K 


2,2,1 


/ i/ \ 
R — 1 (j a I 

£j — [k (uj — a)) 


V w ; 


:) 

£/ = (ac a; J 


2,1,0 


/ p v \ 

\o w y 

Z7 1 

= ft 




I) 

E — a, \iv — 


2,1,1 


/ i/ \ 
it! = 1 a a I 

= uj — a 




*=(s :) 

E — K 


2,0,0 


y a U y 
not existing 


_E not existing 


£/ not 
existing 


1,2,1 


= (a; a) 


*=(::)•*=(: o) 


E = (a 0) 


1,1,0 


(,\ 

it = 1 z/ , .c/ = a 
V a / 




V 17 / 

E — o 




/ \ 

R=\ v , E = uo 




/ ,, \ 






£ = 


1,0,0 


/,\ 

R= a 

\ / 
E not existing 


_E not existing 


-C) 

£ not 
existing 


0,1,0 


i? not existing, 

E = a 


i? not ^ = / « \ 
existing, \ a J 


E = a, R not 
existing 
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Table 3. Submatrices R and E of matrices (5^ 



m', n', A' 


to, n, A 


O 1 1 

2,1,1 


onn 

2,0,0 


1 O 1 

1,Z,1 


9 1 1 


R= {« o) 

E = a, jiu = 


not existing 


D / ,, , A 

r£ — yji V) 

*-(:) 


onn 


/ , , \ 

*=(" S) 


/ ,, ,/ \ 

*=(;: o)' £ - 


K = 1 /i 2/ 1 

£ not 




£< not existing 


existing, /iz/ = 


existing 


1,2,1 










-C/ — ^(X Q; J 


E not existing 


/iz/ — 0; R — a 


1,1,0 


ft — ( L L \ F — rr 

-it — 1 \ . Hi — O 
\ V ) 


R = fi, E not 
existing 


R = ji 

f u \ 


1,1,1 


«-(;)•'•♦ 


R = fi, E not 
existing 


*="■*=( a) 


1,0,0 


«-C) 

E not existing 


R = ji, E not 
existing 


R = fi, 

E not existing 


0,1,0 


R not existing, 

E — a 


i? and E 
not existing 


R not existing 

E = n 



Table 4. Submatrices R and E of matrices /?4 



m', n', A' 


m, n, A 


1,1,0 


1,1,1 


1,0,0 


0,1,0 


1,1,0 


R = H 
E = v 


R = fi 
E = v 


R = fi,E not 
existing 


E — /i,R not 
existing 


1,1,1 


R = H, 
E = v 


R = H 
E = 


R = fi,E not 
existing 


E = fi,R not 
existing 


1,0,0 


R — /i,E not 
existing 


R = fi,E not 
existing 


R — /i,E not 
existing 


i? and E 
not existing 


0,1,0 


E = fi,R not 
existing 


E = fi,R not 
existing 


i? and £ 
not existing 


E = fi,R not 
existing 



These matrices are simplified using the equivalence transformations 



R -»• C/M/ f , £ -> V^^yt. (Al) 

Here [/ and V are unitary matrices whose dimensions are the same as dimensions 
of matrices R and E correspondingly, which satisfy the following relations: UA — 
AU, UB = BV, VC = CU. Transformations (Al) keep equations (21) invariant. 
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B More on the Galilean Rarita-Schwinger equa- 
tion 

Let us prove that equation (55) is consistent and describes a particle with spin s = 3/2. 

Reducing (55) by p m we obtain Am 2 ^ 4 = 0, i.e., if 4 = 0. Whereas, reducing (55) 
by 7m we obtain 

pj n = %p n lJ m - (Bl) 

Finally, comparing (Bl) with equation (55) for index m = 4 we find the following 
consequences of equation (55): 

7nP n ^ = 0, a = 0,1, 2, 3, (B2) 
m^°-p a ^ a = 0, a = 1,2, 3, (B3) 
7o^° + 7 a ^ a = 0, and vj> 4 = 0. (B4) 

On the other hand equation (55) follows from (B2)-(B4), so that equations (55) and 
(B2)-(B4) are equivalent. 

In accordance with (B2) any component of \l/ m satisfies the Levy-Leblond equation 
(compare with Section 2.1). Let us prove now that equations (B2)-(B4) describe 
indeed a particle with spin s=3/2. 

It follows from (B2)-(B4) that, in the rest frame, ^° = ^ 4 = and * a has only 
two non-zero spinor components a = 1,2. Using 7-matrices in realization (14) 
and equation (B4), we conclude that 4f a satisfies the equation 

v a ^ a = 0, • (B5) 

Consequently this function satisfies conditions (31) and (32) as well with s = 3/2. 
This follows from the fact the total spin operator S is a sum of operators of spin one 
and of spin one-half: S a — s a + ^<J a , so that that 

S 2 = — + s • a. (B6) 
4 

Let ^ denotes the column \& 2 , ^ 3 ). In accordance with (B6) the condition 
S 2 \j> = s(s + 1)^ reduces to the form 

^ a - ^e abc a b i> c = 0, (B7) 

for s = 3/2 provided we use the representation with (s a ) bc = itabc, where e abc is a 
totally antisymmetric unit tensor. 

Comparing (B5) with (B7) we conclude that these equations are completely equiv- 
alent since multiplying (B5) by a a we obtain (B7) and multiplying (B7) by s a and 
contracting it with respect to index a we get (B5). 

Thus indeed equation (55) describes a Galilean particle with spin s = 3/2. 



29 



C Contraction of relativistic Proca equation 



Let us consider the relativistic Proca equation (2) and contract it directly to a non- 
relativistic (i.e., Galilei- invariant) approximation. Solutions of this equation are a 
four-vector ^ and a skew-symmetric tensor which transform according to the 
representation -D(|, |) © D(l, 0) © -D(0, 1) of the Lorentz group. 

It was shown in papers [5] and [6] how this representation can be contracted to 
representation -D(3, 1, 1) of the homogeneous Galilei group. This contraction can be 
used to reduce equation (2) to a Galilei-invariant form. To do this it is necessary: 

• To choose the following new dependent variables 

W c = l £ ^c^ bc? B = 

which, in accordance with (2), satisfy the following equations: 
2{p° -K,)R a + p a B + e ahc p b W c = 0, 

(p° + K)N a - e abc p b W c + p a B = 0, ( . 

e abc p b (R c + ±N C ) = KW a , [ } 

\p a N a - p a R a = kB; 

• To act on variables R a , N a , W a and B by a diagonal contraction matrix. This 
action yields the change: 

R a = R a ) N a = e 2 N a , W a = eW a , B = eB, 

where e is a small parameter associated with the inverse speed of light; 

• To change relativistic four-momentum p M and mass k by their Galilean coun- 
terparts p a ,p° and m, where 

p a = e~ 1 p a , p° = po — k and m — -{po + K)e~ 2 ; 

2 

• In each equation in (CI) keep only terms which are multiplied by the lowest 
powers of e. 

As a result we obtain the following equations for R a , N a , W a and B: 
2p°R a + p a B + e abc p b W c = 0, 

e abc p b R c = mW a , (C2) 
p a R a + mB = 0, 

and 

2mN a = e abc p b W c - p a B. (C3) 



30 



The system (C2) is nothing else but the Galilei- invariant equation (1) with ma- 
trices (42) written componentwise. Relation (C3) expresses the extra component N a 
via derivatives of the essential ones, i.e. of W c and B. 

Thus the Galilean analogue of the Proca equation (48) cannot be obtained as a 
non-relativistic limit of the relativistic Proca equation (2) but is a specific modification 
of it. The relativistic counterpart of equation (48) is a specific generalization of (2) 
which will be studied in a separate publication. 
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